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Abstract 

We consider Higgs boson production by gluon fusion in hadron collisions. We 
study the doubly-differential transverse-momentum (q?) and rapidity (y) distribution 
of the Higgs boson in perturbative QCD. In the region of small qx {qr "C Mh, 
Mn being the mass of the Higgs boson), we include the effect of logarithmically- 
enhanced contributions due to multiparton radiation to all perturbative orders. We 
use the impact parameter and double Mellin moments to implement and factorize 
the multiparton kinematics constraint of transverse- and longitudinal-momentum 
conservation. The logarithmic terms are then systematically resummed in exponential 
form. At small qx, we perform the all-order resummation of large logarithms up to 
next-to-next-to- leading logarithmic accuracy, while at large qx (qr ~ Mh), we apply 
a matching procedure that recovers the fixed-order perturbation theory up to next- 
to-leading order. We present quantitative results for the differential cross section in 
qx and y at the LHC, and we comment on the comparison with the qx cross section 
integrated over y. 



This paper is dedicated to the memory of 
Jiro Kodaira, great friend and distinguished colleague 



May 2007 



1 Introduction 



The search for the Higgs boson pQ and the study of its properties (mass, couplings, decay widths) 
at hadron colliders require a detailed understanding of its production mechanisms. This demands 
reliable computations of related quantities, such as production cross sections and the associated 
distributions in rapidity and transverse momentum. In this paper we consider the production of 
the Standard Model (SM) Higgs boson by the gluon fusion mechanism. 

The gluon fusion process gg — > H, through a heavy-quark (mainly, top-quark) loop, is the main 
production mechanism of the SM Higgs boson H at hadron colliders. When combined with the 
decay channels H — > 77 and H — > ZZ, this production mechanism is one of the most important 
for Higgs boson searches and studies over the entire range, 100 GeV ^ M H < 1 TeV, of Higgs boson 
mass Mh to be investigated at the LHC [2] . In the mass range 140 GeV ^ Mh ^ 180 GeV, the gluon 
fusion process, followed by the decay H — > WW — > t+t'vv, can be exploited as main discovery 
channel at the LHC and also at the Tevatron [3], provided the background from tt production 
is suppressed by applying a veto cut on the transverse momenta of the jets accompanying the 
final-state leptons. 

The dynamics of the gluon fusion mechanism in controlled by strong interactions. Detailed 
studies of the effect of QCD radiative corrections are thus necessary to obtain accurate theoretical 
predictions. 

In QCD perturbation theory, the leading order (LO) contribution to the total cross section 
for Higgs boson production by gluon fusion is proportional to ag, as being the QCD coupling. 
The QCD radiative corrections to the total cross section are known at the next-to-leading order 
(NLO) [I]-[7] and at the next-to-next-to-leading order (NNLO) [B]-[I2]- The Higgs boson rapidity 
distribution is also known at the NLO [13] and at the NNLO [T4"l IT5] . The effects of a jet veto 
have been studied up to the NNLO pAj [HJ [15]. We recall that all the results at NNLO have 
been obtained by using the large-M 4 approximation, M t being the mass of the top quark. This 
approximation is justified by the fact that the bulk of the QCD radiative corrections to the total 
cross section is due to virtual and soft-gluon contributions JTHl El EH EH ELT]. The soft-gluon 
dominance also implies that higher-order perturbative contributions can reliably be estimated by 
applying resummation methods [9j of threshold logarithms, a type of logarithmically-enhanced 
terms due to multiple soft-gluon emission. In Ref. [17], the NNLO calculation of the total cross 
section is supplemented with threshold resummation at the next-to-next-to-leading logarithmic 
(NNLL) level; the residual perturbative uncertainty at the LHC is estimated to be at the level 
of better than ±10%. The NNLL+NNLO results [17] are nicely confirmed by the more recent 
computation [H]-[20] of the soft-gluon terms at N 3 LO; the quantitative effect [18] of the additional 
(i.e., beyond the NNLL order) single- logarithmic term at N 3 LO is consistent with the estimated 
uncertainty at NNLL+NNLO. The effect of threshold logarithms on the rapidity distribution of 
the Higgs boson has been considered in Ref. [21]. 

The gluon fusion mechanism at 0(a$) produces a Higgs boson with a vanishing transverse 
momentum qr- A large (or, however, non-vanishing) value of qr can be obtained only starting 
from (9(a|), when the Higgs boson is accompanied by at least one recoiling parton in the final 
state. This mismatch by a power of as is a preliminary indication of the fact that the small-gr 
and large-gr regions are controlled by different dynamics regimes. 
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The large-g^ region is identified by the condition qx ~ M#. In this region, the perturbative 
series is controlled by a small expansion parameter, as(Mfj), and calculations based on the trunca- 
tion of the series at a fixed order in as are theoretically justified. The LO, i.e. C?(a|), calculation 
is reported in Ref. [22]. The results of Ref. [22] and the higher-order studies of Refs. (23], [24] show 
that the large- M t approximation is sufficiently accurate also in the case of the qx distribution when 
qr^Mn, provided qx^M t . Using the large- M t approximation, the NLO QCD computation of 
the qx distribution of the SM Higgs boson is presented in Refs. [251 I2S1 [271 HH EH]- QCD correc- 
tions beyond the NLO are evaluated in Ref. [28], by implementing threshold resummation at the 
next-to-leading logarithmic (NLL) level. The results of the numerical programs of Refs. [TH [15] 
can also be safely (i.e. without encountering infrared divergences) extended from large values of 
qT to qx = 0: in the small-g^ region these programs evaluate the qx distribution up to NNLO. 

In the small-g^ region [qx <C M#), where the bulk of events is produced, the convergence 
of the fixed-order expansion is definitely spoiled, since the coefficients of the perturbative series 
in as(Mfi) are enhanced by powers of large logarithmic terms, \n m (Mfj / q T ) . The logarithmic 
terms are produced by multiple emission of soft and collinear partons (i.e. partons with low trans- 
verse momentum). To obtain reliable perturbative predictions, these terms have to be resummed 
to all orders in as- The method to systematically perform all-order resummation of classes of 
logarithmically-enhanced terms at small qx is known |29]-[36]. In the case of the SM Higgs boson, 
resummation has been explicitly worked out at leading logarithmic (LL), NLL [35l [37] and NNLL 
[38] level. 

The fixed-order and resummed approaches at small and large values of qx can then be matched 
at intermediate values of qx, to obtain QCD predictions for the entire range of transverse momenta. 
Phenomenological studies of the SM Higgs boson qx distribution at the LHC have been performed 
in Refs. [39]-|45j. by combining resummed and fixed-order perturbation theory at different levels 
of theoretical accuracy. Other recent studies of various kinematical distributions of the SM Higgs 
boson at the LHC are presented in Refs. [4"6]-[50j. 

In Refs. [HI S3] we studied the Higgs boson q T distribution integrated over the rapidity. In the 
small-g^ region, the logarithmic terms were systematically resummed in exponential form by work- 
ing in impact-parameter and Mellin-moment space. A constraint of perturbative unitarity was 
imposed on the resummed terms, to the purpose of reducing the effect of unjustified higher-order 
contributions at large values of qr and, especially, at intermediate values of qx- This constraint 
thus decreases the uncertainty in the matching procedure of the resummed and fixed-order con- 
tributions. Our best theoretical predictions were obtained by matching NNLL resummation at 
small qx and NLO perturbation theory at large qx- NNLL resummation includes the complete 
NNLO result at small qr, and the unitarity constraint assures that the total cross section at 
NNLO is recovered upon integration over qx of the transverse-momentum spectrum. Considering 
SM Higgs boson production at the LHC, we concluded [33] that the residual perturbative QCD 
uncertainty of the NNLL+NLO result is uniformly of about ±10% from small to intermediate 
values of transverse momenta. 

In this paper we extend our study to include the dependence on the rapidity of the Higgs 
boson. Using the impact parameter and double Mellin moments, we can perform the extension 
by maintaining all the main features of the resummation formalism of Refs. [361 [33]- We are then 
able to present results up to NNLL+NLO accuracy for the doubly-differential cross section in qx 
and rapidity at the LHC. 
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The paper is organized as follows. In Sect. [2] we recall the main aspects of the resummation 
formalism, and we illustrate the steps that are necessary to include the dependence on the rapidity 
in the q F resummed formulae. In Sect. [3] we apply the formalism to the production of the SM 
Higgs boson at the LHC, and we perform quantitative studies on the qx and rapidity dependence 
of the doubly-differential cross section. Some concluding remarks are presented in Sect. HI Addi- 
tional technical details on the double Mellin moments of the resummation formulae are given in 
Appendix lAl 



2 Rapidity dependence in qx resummation 

We consider the inclusive hard-scattering process 

hM + h 2 (p 2 ) ^ H(y,q T ,M H ) + X , (1) 

where the collision of the two hadrons hi and h 2 with momenta p\ and p 2 produces the Higgs 
boson H, accompanied by an arbitrary and undetected final state X. The centre-of-mass energy 
of the colliding hadrons is denoted by The rapidity, y, of the Higgs boson is defined in the 
centre-of-mass frame of the colliding hadrons, and the forward direction {y > 0) is identified by 
the direction of the momentum p\. 

According to the QCD factorization theorem, the doubly-differential cross section for this 
process is 



da 



dy dq 



T (y,q T ,M H ,s) = dxi \ dx<2 /ai/fci( x i>A4) fa 2 /h 2 (x 2 ,n 2 F ) 

T „. „„ JO JO 

(y,q T ,M H ,s\a s (v 2 R ),fJ? R ,fi 2 F ) , (2) 



11,12 

da„ 



dy dq T 

where f a /h(x,fi F ) (a = qf,q~f,g) are the parton densities of the colliding hadrons at the factoriza- 
tion scale fi F , da a b are the partonic cross sections, and /i# is the renormalization scale. Throughout 
the paper we use parton densities as defined in the MS factorization scheme, and ces(q 2 ) is the 
QCD running coupling in the MS renormalization scheme. The rapidity, y, and the centre-of-mass 
energy, s, of the partonic cross section (subprocess) are related to the corresponding hadronic vari- 
ables y and s: 

V = V - 2 ' S = XlX2S ' ( 3 ) 
with the kinematical boundary \y\ < In ^Js/M 2 ( \y\ < In a/s/M 2 ) and s > M 2 (s > M 2 ). 

The partonic cross section da a b is computable in QCD perturbation theory. Its power series 
expansion in as contains the logarithmically-enhanced terms, (ag/gfO ln m (M| r /gy), that we want 
to resum. To this purpose, we use the general (process-independent) strategy and the formalism 
described in detail in Ref . [H] . The only difference with respect to Ref . [44J is that the resummation 
is now performed at fixed values of the rapidity y, rather than after integration over the rapidity 
phase space. In the following we briefly recall the main steps of the resummation formalism, and 
we point out explicitly the differences with respect to Ref. [44] . 
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We first rewrite (see Sect. 2.1 in Ref. [2]) the partonic cross section as the sum of two terms, 



j~( res -) j"(fi n -) 
dy dq\ dy dq\ dy dq\ 



+ (4) 



The logarithmically-enhanced contributions are embodied in the 'resummed' component dai^- 
The 'finite' component dd^al is free of such contributions, and it can be computed by truncation 
of the perturbative series at a given fixed order (LO, NLO and so forth). In practice, after having 
evaluated da aia2 and its resummed component at a given perturbative order, the finite component 
da^a} is obtained by the matching procedure described in Sects. 2.1 and 2.4 of Ref. 



The resummation procedure of the logarithmic terms has to be carried out [30]-[34J in the 
impact-parameter space, to correctly take into account the kinematics constraint of transverse- 
momentum conservation. The resummed component of the partonic cross section is then obtained 
by performing the inverse Fourier (Bessel) transformation with respect to the impact parameter 
b. We writeS 



M f°° b 

(y,q T ,M H ,s;a s ) = — ^ / db - J (bq T ) W ai a 2 (y,b, M H , s; a s ) , (5) 

Liy LLq T S Jq Z 

where Jq(x) is the Oth-order Bessel function, and the factor W embodies the all-order dependence 
on the large logarithms ln(M#5) 2 at large b, which correspond to the g^-space terms \n(M H /q^) 
(the limit qx <C Mh corresponds to Mnb ^> 1, since b is the variable conjugate to qx). 

In the case of the gr cross section integrated over the rapidity, the resummation of the large 
logarithms is better expressed [361 14*4"] by defining the iV-moments Wat of W with respect to 
z = M H /s at fixed Mh- In the present case, where the rapidity is fixed, it is convenient (see 
e.g. Refs. [511 [52]) to consider 'double' (JVi, A^-moments with respect to the two variables z\ = 
e + yM H /V§ and z 2 = e^M^/VI at fixed M H (note that < z { < 1). We thus introduce W {Nl ' N2) 
as follows: 



WZf*\b,M H ;a s ) = f'dz.z^- 1 f dz 2 z^~ 1 W aia2 (y,b, M H ,s;a s ) 

Jo Jo 



(6) 



More generally, any function h(y;z) of the variables y (\y\ < — hx^/z) and z (0 < z < 1) can 
be considered as a function of the two variables z\ = e +y ^fz and z 2 = e~ y y/z. Thus, throughout 
the paper, the (Ni, A^ 2 )-moments h^ Nl ' N2 ^ of the function h(y; z) are defined as 

h ( Nl ,N 2 ) _ I ^ z Nr-l I ^ Z N 2 -1 h ( y . z) ^ where . y = 1 ln £l ^ = Zi ^ _ (7) 

Jo Jo 2 z 2 

Note that the double Mellin moments can also be obtained (see e.g. Ref. [53]) by introducing 
a Fourier transformation with respect to y (with conjugate variable v = i(N 2 — N±)) and then 
performing a Mellin transformation with respect to z (with conjugate variable iV = (N± + N 2 )/2): 

h (N u N 2 ) = / dz z n-i I dy e iu V h ^ y . ^ ^ where :Nl = N + iv j 2 ,N 2 = N - iv/2 . (8) 

Jo J-oo 

tin the following equations, the functional dependence on the scales fiR and up is understood. 
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The convolution structure of the QCD factorization formula ([2]) is readily diagonalized by 
considering (Ni, iV^-moments: 

do^iM = f ai/hl ,N 1+ t fa M+ l da^ , (9) 

11,12 

where f a /h,N = J dx x N ~ l f a /h( x ) are the customary iV-moments of the parton distributions. 

The use of Mellin moments also simplifies the resummation structure of the logarithmic terms 
in ddaia'2 ( - Nl,N2 - ) . The perturbative factor VV^^'^ can indeed be organized in exponential form as 
follows: 

W^\b,M H ;a s ) = H^\M H ,a s ) exp{0W^(a s , L)} , (10) 

where 

b = 2e 7B (■je = 0.5772 ... is the Euler number) and, to simplify the notation, the dependence 
on the flavour indeces has been understood. 

The structure of Eq. ffTUl) is in close analogy to the cases of soft-gluon resummed calculations 
for hadronic event shapes in hard-scattering processes [M] and for threshold contributions to 
hadronic cross sections [55l |5T| [56]. The function T-c( Nl < N2 ) (which is process dependent) does not 
depend on the impact parameter b and, therefore, its evaluation does not require resummation of 
large logarithmic terms. It can be expanded in powers of as as 

H^'^ (M H , as) = a (a s ,M H ) [l + ^7^*0(1) + f^V ftC^.iW) + _ 1 ? (12) 

L 7T V 7T / J 

where <Jo(as, Mh) is the lowest-order partonic cross section for Higgs boson production. The 
form factor exp{(?} is process independent^, it includes the complete dependence on b and, in 
particular, it contains all the terms that order-by-order in as are logarithmically divergent when 
b — > oo. The functional dependence on b is expressed through the large logarithmic terms a§L m 
with 1 < m < In. More importantly, all the logarithmic contributions to Q with n + 2 < m < 2n 
are vanishing. Thus, the exponent Q can systematically be expanded in powers of as, at fixed 
value of A = asL, as follows: 

g^\as,L) = Lg^\a s L)+g^^\a s L) + ^ g® ^\a s L) + . . . . (13) 

7T 

The term Lg^ collects the leading logarithmic (LL) contributions «gl/ n+1 ; the function g^ re- 
sums the next-to-leading logarithmic (NLL) contributions a^L n ; g( 3 ' controls the next-to-next-to- 
leading logarithmic (NNLL) terms agL n_1 , and so forth. 

Note that we use the logarithmic variable L (see Eq. (fTTj) ) to parametrize and organize the 
resummation of the large logarithms ln(M#6) 2 . We recall the main motivations [H] for this choice. 
In the resummation region Mub 3> 1, we have L ~ ln(M#6) 2 and the use of the variable L is 
fully legitimate to arbitrary logarithmic accuracy. When M^b < 1, we have L — > (whereat 



■fMore precisely, it depends only on the flavour of the colliding partons (see Appendix |A|) . 

§ As shown in Appendix B of Ref. [44] (see Eqs. (131) and (132) therein), after inverse Fourier transformation 
to qr space, the 6-dependent functions ln™(Mjy&) 2 and L n lead to quite different behaviours at large qr- When 
qr ^> Mh, the behaviour (1/q 2 -) hi n ~ l {qT /Mr) (which is not integrable when qr — ► oo) produced by ln"(M#6) 2 
is damped (and made integrable) by the extra factor <J qx /Mh exp ( — 6 q? / Mh ) produced in the case of L n . 
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\ia(M H b) 2 — > oo !) and exp{(?(a s , L)} — > 1. Therefore, the use of L reduces the effect produced 
by the resummed contributions in the small-6 region (i.e., at large and intermediate values of 
qx), where the large-6 resummation approach is not justified. In particular, setting 6 = (which 
corresponds to integrate over the entire qx range) we have exp{(?(as, L)} = 1: this property 
can be interpreted [S] as a constraint of perturbative unitarity on the total cross section; the 
dynamics of the all-order recoil effects, which are resummed in the form factor exp{Q(as, L)}, 
produces a smearing of the fixed-order q? distribution of the Higgs boson without affecting its 
total production rate. 

The resummation formulae (jlOp . fTl2|) and ffl3|) can be worked out at any given (and arbitrary) 
logarithmic accuracy since the functions 7i and Q can explicitly be expressed (see Ref. [H]) in 
terms of few perturbatively-computable coefficients denoted by A^ n \ B^ n \ H^ n \ Cjy , jffi. The key 
role of these coefficients to fully determine the structure of transverse-momentum resummation was 
first formalized by Collins, Soper and Sterman [3U [32], [36] . The present status of the calculation 
of these coefficients for Higgs boson production is recalled in Sect. [3j 

In the case of the cross section integrated over the rapidity, Eq. (flUj) is still valid, provided the 
double (Nx, N 2 ) -moments are replaced by the corresponding single iV-moments Wn,'Hn,Gn (see 
Sect. 2.2 in Ref. |44|). The relation between double and single moments can easily be understood 
by inspection of Eqs. (jSJ)-©- We see that setting v = in Eq. (jSJ) is exactly equivalent to integrate 
the cross section over the rapidity. Therefore, the functions Wn, T~Cn, Gn in Ref- pH] are obtained 
by simply setting Nx — N% — N in the corresponding functions VV^ 1 '^, ri^ Nl ' N2 \ Q( Nl > N ti of 

Eq. (unj. 

Moreover, from the results presented in Ref. [H] , we can straightforwardly obtain the functions 
q-^{N 1 ,N 2 ) anc j g{Ni,N 2 ) f rom fa e functions ri^ and Gn- Roughly speaking, we simply have 

G {N ^ ] = \ (Gm + Gn 2 ) , rt N ^ = [H Nl H N2 ] 1/2 . (14) 

More precisely, these equalities are valid in the simplified case where there is a single species of 
partons (e.g. only gluons). In the following we comment on the physical picture that leads to 
Eq. (I14p . The generalization to considering more species of partons does not require any further 
conceptual steps: it just involves algebraic complications related to the treatment of the flavour 
indeces. The multiflavour case is briefly illustrated in Appendix [A3 

In the small-g^ (large-6) region that we are considering, the kinematics of the Higgs boson is 
fully determined by the radiation of soft and collinear partons from the colliding partons (hadrons) 
in the initial state. The radiation of soft partons cannot affect the rapidity of the Higgs bosons. On 
the contrary, the radiation of partons that are collinear to pi (j>2), i.e. in the forward (backward) 
region, decreases (increases) the rapidity of the Higgs boson as a consequence of longitudinal- 
momentum conservation (see Eq. ([3])). Since the emissions of collinear partons from px and p 2 
are dynamically uncorrelated (factorized from each other), correlations arise only from kinemat- 
ics. The use of the (Nx, AT 2 ) -moments exactly factorizes (see Eqs. (T5]) and (Q) the kinematical 
constraint of longitudinal-momentum conservation. It follows that the (Nx, N2) -dependence of 
W(iVi,iV 2 ) ig given by the proc i U ct of two functions (say, W {Nl ' N2) = M[ Nl) M 2 N2) ) that depends 
only on Nx or N 2 , respectively. If all the partons have the same flavour, the two functions should 
be equal, and Eq. (HU) directly follows from [W^ Nl!N ^] Nl=N2=N = W N . 

The formalism illustrated in this section defines a systematic 'order-by-order' (in extended 
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sense) expansion [H] of Eq. (HI) : it can be used to obtain predictions with uniform perturbative 
accuracy from the small-g^ region to the large-g^ region. The various orders of this expansion 
are denoted^ as LL, NLL+LO, NNLL+NLO, etc., where the first label (LL, NLL, NNLL, . . . ) 
refers to the logarithmic accuracy at small qx and the second label (LO, NLO, . . . ) refers to 
the customary perturbative ordeiHH at large qx- To be precise, the NLL+LO term of Eq. (J4j) is 
obtained by including the functions and the coefficient TiS 1 ' (see Eqs. ( TTBl and ( |T2l) ) in 

the resummed component, and by expanding the finite (i.e. large-gj 1 ) component up to its LO 
term. At NNLL+NLO accuracy, the resummed component includes also the function g N and the 
coefficient H {2) (see Eqs. ([TBI and ffl2l) ). while the finite component is expanded up to NLO. It 
is worthwhile noticing that the NNLL+NLO (NLL+LO) result includes the full NNLO (NLO) 
perturbative contribution in the small-gy region. 

We recall [S] that, due to our actual definition of the logarithmic parameter L in Eq. (flOl) and 
to our matching procedure with the perturbative expansion at large qx, the integral over qx of the 
qx cross section exactly reproduces the customary fixed-order calculation of the total cross section. 
This feature is not affected by keeping the rapidity fixed. Therefore, the NNLO (NLO) result for 
total cross section at fixed y is exactly recovered upon integration over qx of the NNLL+NLO 
(NLL+LO) qx spectrum at fixed y. 

Within our formalism, resummation is directly implemented, at fixed Mh, in the space of 
the conjugate variables N\,N2 and b. To obtain the cross section in Eq. (T5]), as function of the 
kinematical variables s, y and qx, we have to perform inverse integral transformations. These 
integrals are carried out numerically. We recall jH] that the resummed form factor (i.e., each of 
the functions g^ k \a^L) in Eq. (fl3|) ) is singular at the value of b where asifJ^L = n/p (/3 is the 
first-order coefficient of the QCD (3 function). This singularity has its origin from the presence of 
the Landau pole in the running of the QCD coupling a^{q 2 ) at low scales. When performing the 
inverse Fourier (Bessel) transformation with respect to the impact parameter b (see Eq. ((5])), we 
deal with this singularity by using a 'minimal prescription' [561 [57]: the singularity is avoided by 
deforming the integration contour in the complex b space (see Ref. [5 TJ ) . We note that the position 
of the singularity is completely independent of the values of Aq and A^. Thus, the inversion of 
the Mellin moments is performed in the customary way (in Mellin space there are no singularities 
for sufficiently- large values of Re Aq and Re A 2 ). In this respect, going from single A-moments 
(as in Ref. [H]) to double (Aq, iV^-moments (as in the present case, where the rapidity is kept 
fixed) is completely straightforward, with no additional (practical or conceptual) complications. 



3 Higgs boson production at the LHC 

In this section we apply the resummation formalism of Sect. [2] to the production of the Standard 
Model Higgs boson at the LHC. We closely follow our previous study of the single differential 
(with respect to qx) cross section, with the same choice of parameters as stated in Sect. 3 of 

^In the literature on qx resummation, other authors sometime use the same labels (NLL, NLO and so forth) 
with a meaning that is different from ours. 

"We recall that the LO term at small qr (i.e. including the region where qx — 0) is proportional to «g, whereas 
the LO term at large qx is proportional to ag . This mismatch of one power of as (and the ensuing mismatch of 
notation) persists at each higher order (NLO, NNLO, ...). 
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Ref. [44]. Therefore, the integration over y of the double differential (with respect to y and q?) 
cross sections presented in this section returns the qx cross sections of Ref. [44]. As a cross-check 
of the actual implementation of the calculation, we have verified that after integration over the 
rapidity the numerical results in Ref. [44] are reobtained within a high accuracy. 

As in Refs. [TT1I44] . we use an 'improved version' [16] of the large- M t approximation. The cross 
section is first computed by using the large- M t approximation. Then, it is rescaled by a Born level 
factor, such as to include the exact lowest-order dependence on the masses, M t and M&, of the top 
and bottom!**! quarks, which circulates in the heavy-quark loop that couples to the Higgs boson. 
We use the values M t = 175 GeV and M b = 4.75 GeV. As discussed in Ref. [17] and recalled in 
Sect. [[J this version of the large- M t approximation is expected to produce an uncertainty that is 
smaller than the uncertainties from yet uncalculated perturbative terms from higher orders. 

For the sake of brevity, we present quantitative results only at NNLL+NLO accuracy, which 
is the highest accuracy that can be achieved by using the present knowledge of exact perturbative 
QCD contributions (resummation coefficients and fixed-order calculations [25]- [2?]). We use the 
MRST2004 set [58] of parton distribution functions at NNLO. The use of NNLO parton densities 
consistently matches the NNLL (NNLO) accuracy of our partonic cross section in the region of 
small and intermediate values of qx- 

Resummation up to the NLL level is under control from the knowledge of the perturbative 
coefficients [35] and H {1) [37J. To reach the NNLL+NLO accuracy, the form factor 

function g( Nl > N2 > in Eq. (TIB"]) must include the contribution from ( Nl > N2 ^ (which is controlled 
by the coefficients B (2) [38] and A® [59]), and the coefficient function H {Nl > N2) in Eq. <$2ft has 
to be evaluated up to its second-order term TiP^ ( Nl > N2 \ In Ref. [44] we exploited the unitarity 
constraint Q(a§, L)\b=o = to numerically derive an approximated form of the coefficient ri^ 
from the NNLO calculation [12] of the total cross section. The recent calculation of Ref. |15j . 
which is based on the complete evaluation of US 2 ' ( Nl > N2 > in analytic form, allows us to gauge the 
quality of the approximated form. We find that the use of the ri^ of Ref. [44] leads to differences 
of about 1% with respect to the exact computation of the rapidity cross section at NNLO. 

All the numerical results in this section are obtained by fixing the renormalization and fac- 
torization scales at the value = /ip = Mjj. The 'resummation scale' Q (the auxiliary scale 
introduced in Ref. [44] to gauge the effect of yet uncalculated logarithmic terms at higher orders) 
is also fixed at the value Q = Mh- The mass of the Higgs boson is set at the value Mjj = 125 GeV. 

We start our presentation of the predictions for Higgs boson production at the LHC by consid- 
ering the qx dependence of the cross section at fixed values of the rapidity. In Fig. [TJ we set y — 
and we compare the customary (when q? > 0) NLO calculation (dashed line) with the resummed 
NNLL+NLO calculation (solid line). 

As expected, the NLO result diverges to — oo as qx — > and, at small values of qx, it has an 
unphysical peak that is produced by the numerical compensation of negative leading logarithmic 
and positive subleading logarithmic contributions. The presence of this peak is not accidental. At 

**We note that the Born level cross section is not insensitive to the contribution of the bottom quark. Adding 
the bottom-quark loop to the top-quark loop in the scattering amplitude produces a non-negligible interference 
effect in the squared amplitude. The relative effect of the bottom quark decreases the Born level cross section by 
about 11% if M H = 125 GeV, and by about 3% if M H = 300 GeV. If M H >500 GeV, the relative effect of the 
bottom quark is always smaller than 1%. 
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M H =125 GeV y=0 




20 40 60 80 100 



q T (GeV) 

Figure 1: The qr spectrum at the LHC with Mh = 125 GeV and y = 0: results at NNLL+NLO 
(solid line) and NLO (dashed line) accuracy. The inset plot shows the ratio K (see Eq. 173]) ,) of 
the corresponding qx cross sections, fixing y = (solid line) and integrating them over the full 
rapidity range (dashed line). 

large qT, the perturbative expansion at any fixed order has no pathological behaviour: it leads to 
a positive cross section, whose value decreases as qT increases. When qT — > 0, instead, any fixed- 
order calculation diverges alternatively to ±oo depending on the perturbative order. Therefore, to 
go smoothly from the large-g^ behaviour to the small-g^ limit, the NLO (or N 3 LO, and so forth) 
calculation of the cross section has to show at least one peak in the intermediate-gT region. 

We recall once more that the label NLO in Fig. prefers to (and originates from) the perturbative 
expansion at large qT- To avoid possible misunderstandings (coming from such a label) when 
interpreting the dashed (NLO) curve in the small-gr region, we point out that, the only difference 
produced in Fig. [I] by the NNLO calculation at small qT (this calculation can be carried out, for 
example, by using the NNLO codes of Refs. [HI [15]) is a spike around the point qT = 0. More 
precisely, as long as qT ^ 0, the dashed curve is exactly the result of the NNLO calculation of the 
qT cross section at small qT- The only difference introduced in the plot by this NNLO calculation 
would occur in the first bin (with arbitrarily small size) that includes the point qT = 0. The 
NNLO value of the qT cross section in this first bin is positive and fixed by the value of the NNLO 
total cross section^. Of course, owing to the increasingly negative behaviour of the qT distribution 
when q T — > 0, the NNLO value of the qT cross section in the first bin increases by decreasing the 
size of that bin. 

The resummed NNLL+NLO result in Fig. [T] is physically well-behaved at small q? (it vanishes 
t^By definition, the integral over qr of d 2 a/ (dqx dy) at NNLO is equal to da/dy at NNLO. 



9 



as qx — > and has a kinematical peak at ~ 12 GeV), and it converges to the expected NLO 
result only when is definitely large (gr — Mff)- 

To quantify more clearly the effect of the resummation on the the NLO result, the value at 
y = of the qx dependent K-factor, 



K(q T ,y) 



da 



NNLL+NLO 



I (dq T dy) 



(15) 



da NLO /(dq T dy) 

is shown in the inset plot of Fig. [TJ The dashed line shows the analogous K-factor as computed 
from the ratio of the rapidity integrated cross sections. The similarity between these two K-factors 
is a first indication of the mild rapidity dependence of the resummation effects. By inspection 
of the inset plot, we note that NNLL resummation is relevant not only at small qx, but also in 
the intermediate-g^ region: as soon as q? < 80 GeV, the resummation effects are larger than 20%. 
Of course, the fact that K ~ 1 at q? ~ 24 GeV is purely accidental: it simply follows from the 
unphysical behaviour of the fixed-order perturbative expansion at small qx- 
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Figure 2: The q T spectrum at the LHC with M H = 125 GeV and y = 2: results at NNLL+NLO 
(solid line) and NLO (dashed line) accuracy. The inset plot shows the ratio K (see Eq. 173]) ,) of 
the corresponding qx cross sections, fixing y = 2 (solid line) and integrating them over the full 
rapidity range (dashed line). 

Considering other values of the rapidity, from the central to the off-central rapidity region, we 
find the same features as observed at y — 0. Our results of the qx spectrum at y — 2 are presented 
in Fig. [2j The NNLL+NLO spectrum has a peak at qr ~ H GeV. As happens in the case of 
the qx distribution integrated over y, the effect of NNLL resummation is definitely non-negligible 
starting from relatively-high values of qT- For example, at qx = 50 GeV the NNLL+NLO result 
is about 30% higher than the NLO result. 
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To analyze the rapidity dependence in more detail, we study the doubly-differential cross 
section at fixed values of qx- In Figs. [3]and|U we show quantitative results at two typical values 
of the transverse momentum, qx = 15 GeV and qx = 40 GeV, in the small-g^ and intermediate-g^ 
region, respectively. 
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Figure 3: The rapidity spectrum at the LHC with M H = 125 GeV and qx = 15 GeV: results at 
NNLL+NLO (solid line) and NLO (dashed line) accuracy. The inset plot shows the K-factor as 
defined in Eq. (fl5]). 

Figure [3] shows the rapidity distribution at NNLL+NLO (solid line) and NLO (dashes) accuracy 
when qx = 15 GeV. At this value of qx, the effect of NNLL resummation reduces the cross section. 
For example, when y = the reduction effect is about 25%. As can be observed in the inset plot, 
the relative contribution from the resummed logarithmic terms is rather constant in the central 
rapidity region, and its dependence on y only appears in forward (and backward) region, where 
the cross section is quite small. 

When qx = 40 GeV (see Fig. H]), instead, the effect of NNLL resummation increases the 
absolute value of the cross section. For example, when y = the NLO cross section is increased 
by about 22%. Nonetheless, as for the relative effect of resummation and the rapidity dependence 
of the K-factor, we observe features that are very similar to those in Fig. [31 The resummation 
effects have a very mild dependence on y in the central and (moderately) off-central regions, and 
this explains the remarkable similarity between the solid and dashed lines in the inset plot of Figs. [1] 
and El Since the kinematical region where \y \ ^2 accounts for most of the total cross section, when 
comparing the ratio K(qx,y) to the analogous ratio of the y-integrated cross sections, hardly any 
differences are expected, unless the large-rapidity region is explored. 

The mild rapidity dependence of the qx shape of the resummed results can be studied with a 
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Figure 4: The rapidity spectrum at the LHC with Mh = 125 GeV and qr = 40 GeV: results at 
NNLL+NLO (solid line) and NLO (dashed line) accuracy. The inset plot shows the K-factor as 
defined in Eq. 173]) . 

finer resolution by denning the following ratio: 

d 2 a/(dq T dy) 

R{qT]y)= da/d qT • (16) 

This ratio gives the doubly-differential cross section normalized to the qx cross section integrated 
over the full rapidity range. For comparison, we consider also the ^-integrated version of the 
cross section ratio in Eq. (jl6p . and we define the ratio 

R y = d -^ (17) 
a 

of the rapidity cross section da/dy over the total cross section a. 



We have computed the ratio in Eq. ( 1161) by using the resummed qr cross sections at NNLL+NLO 
accuracy. The results, as a function of q?, are presented in Fig. [5] (solid lines) at two different 
values, y = and y — 2, of the rapidity. The results of the analogous ((/^-independent) ratio R y 
(computed^] at NNLO with the numerical programs of Refs. [HI [15]) at the corresponding values 
of rapidity are also reported (dotted lines) in Fig. [5J The dashed lines in Fig. [5] correspond to 
the computation of Eq. ffTBl) by using the qx cross sections at NLO: we see that the dashed and 
solid lines are very similar (as expected from the similarity of the dashed and solid lines in the 



ttThe numerical accuracy of this computation is better than about 2%-3%. Owing to the unitarity constraint 
in our resummation formalism, the same result (with a similar numerical accuracy) can be obtained by integration 
over qt of the resummed qx cross sections. 
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Figure 5: The rescaled qx spectrum (as defined by the ratio R(qr] y) in Eq. ICto)) ) at the LHC with 
Mn — 125 GeV. The solid (dashed) lines correspond to the NNLL+NLO (NLO) results at two 
different values of the rapidity: y = (upper) and y = 2 (lower). The dotted lines refer to the 
corresponding values of the ratio R y (see Eq. (T7\)). 

inset plot of Figs. CD and [2]). As discussed below, the results in Fig. [5] show that the cross section 
decreases and the qr spectrum softens when the rapidity increases. 

We observe that the lines at y — lie above the lines at y = 2; this is just a consequence of 
the fact that the cross sections (both at fixed qx and after integration over qx) decrease when y 
increases. 

At fixed y, R(qr] y) is not constant: it depends (though very slightly) on q T . We note that the 
corresponding upper and lower lines in Fig. [5] have different slopes with respect to qx'- fixing q?, 
the qx slope of R(qr] y) decreases from positive to negative values as y increases from y = to 
y = 2, thus showing that the qx spectrum becomes slightly softer at larger rapidity. 

In general, as \y\ increases, the hardness of the qx shape of d 2 o~/(dqT dy) decreases. Since the 
cross section decreases by increasing the rapidity, the hardness of da/dqx (the denominator in 
Eq. ffl6|) ) is intermediate between the values of the hardness of d 2 a/(dqTdy) (the numerator in 
Eq. ffl6l) ) at y = and at large \y\. As a consequence, the qx slope of R(qx',y) is necessarily 
positive when y = 0. Note that the qx slope is already negative when y = 2 (Fig. [5]): this is a 
consequence of the fact that the bulk of the cross section is in the rapidity region \y\ <2. 

Our qualitative illustration of the results in Fig. can be accompanied by some quantitative 
observations. We note that the rapidity dependence of the cross sections is sizeable: going from 
y — to y — 2, the ratio Ry decreases by about 43%; comparable variations affect the ratio 
R(qt] y), which is not very different from R y and it is slowly dependent on q T . Indeed, at fixed y, 
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the ratio R(qx', y) at NNLL+NLO accuracy has a small and nearly constant slope from low values 
of qx around the peak (say, qx ~ 10 GeV) to qx = 100 GeV; varying q T in this region, R(qx] y) 
increases by about 11% when y — 0, and it decreases by about 16% when y = 2. In the same 
range of qx and y, the values of R(qx',y) at NNLL+NLO (solid lines) and at NLO (dashed lines) 
are very similar: although this is expected at large qx, the differences never exceed the level of 
about 4% even at values of qx as low as qx ~ 10 GeV. 

In summary, the results in Fig. [5] show that, when \y\ increases from the central to the (mod- 
erately) off-central region, the cross sections vary more in absolute value than in qx shape. These 
features deserve some words of discussion. 

We first consider the total cross section a and the rapidity cross section da/dy. We recall (see 
Sect. CD) that the value of these cross sections is sizeably affected by QCD radiative corrections. 
The bulk of the effect is due to the radiation of virtual and soft gluons, and they cannot affect the 
rapidity of the Higgs boson. As a consequence, the ratio R y has little sensitivity to perturbative 
QCD corrections. The decreases of R y as \y\ increases is mainly driven by the decrease of the gluon 
density f g (x, Mjj) as x increases. Considering the large-g^ region, similar arguments apply to the 
qx cross sections da/dqx and da/(dqxdy), and similar conclusions apply to the ratio R(qx,y)- In 
the small-gr region, we have to consider the additional and large effect produced on the qx cross 
sections by the logarithmically-enhanced terms \a m {Mjj / q T ) . These terms are due to the radiation 
of soft and collinear partons. As already discussed in Sect. El the rapidity of the Higgs boson can 
be varied only by collinear radiation, while soft radiation can only lead to on overall (independent 
of y) rescaling of the qx cross sections. At the LL level, only soft radiation contributes (the LL 
function g^ x > in Eq. (1131) does not depend on Ni and N 2 ) and all the logarithmic terms cancel in 
the ratio R(qx, y)- The y sensitivity of R(qx, y) starts at the NLL level. The corrections produced 
on the dominant soft-gluon effects by the collinear radiation are physically [29] well approximated 
by varying the scale fi of the gluon density from fi ~ M H to fi ~ qx- As a consequence, the 
variations of the hardness of the of the qx cross sections are mainly driven by dlnf g (x, q T )/d\nq T , 
the amount of scaling violation of the gluon density. Since the scaling violation decreases as x 
increases, the hardness of da/(dq T dy) decreases and the qx spectrum softens as \y\ increases. Note 
that, by increasing x, the gluon density decreases faster than its scaling violation: this explains 
why daj[dqxdy) varies more in absolute value than in qx shape when \y\ increases. 

We conclude this section with some comments about the theoretical uncertainties on the 
doubly-differential cross section da/{dq T dy) at NNLL+NLO accuracy. In Ref. [H] the pertur- 
bative QCD uncertainties on da/dqx were investigated by comparing the results at NNLL+NLO 
and NLL+LO accuracies and by performing scale variations at NNLL+NLO level. We also con- 
sidered the inclusion of non-perturbative contributions, and we found that they lead to small 
corrections provided qx is not very small. From these studies we concluded that the NNLL+NLO 
result has a QCD uncertainty of about ±10% in the region from small (around the peak of the 
qx distribution) to intermediate (say, roughly, g T <M#/3) values of transverse momenta. Similar 
studies can be carried out in the case of the doubly-differential cross section da/(dq T dy). These 
studies are not reported here since their results and the ensuing conclusions are very similar to 
those in Ref. [H] . The reason for this similarity is a feature that we have pointed out throughout 
this section: the qx resummation effects have a very mild dependence on the rapidity and, thus, 
they are almost unchanged when comparing da/{dqxdy) with da/dqx (equivalently, they largely 
cancel in the ratio in the cross section ratio of Eq. ffTBT) ). 
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4 Summary 



We have considered the resummation of the logarithmically-enhanced QCD contributions that 
appear at small transverse momenta when computing the spectrum of a Higgs boson produced 
in hadron collisions. In our previous work on the subject [3H 03], the rapidity of the Higgs boson 
was integrated over: resummation was implemented by using a formalism based on a transform 
to impact parameter and Mellin moment space. In this paper we have extended the resummation 
formalism to the case in which the rapidity is kept fixed, and we have considered the doubly- 
differential cross section with respect to the transverse momentum and the rapidity. We have 
shown that this extension can be carried out without substantial complications: it is sufficient to 
enlarge the conjugate space by introducing a suitably-defined double Mellin transformation. 

The main aspects of our method [36], 133], which are recalled here, are unchanged by the inclu- 
sion of the rapidity dependence. The resummation is performed at the level of the partonic cross 
section, and the parton densities are factorized as in the customary fixed-order calculations. The 
formalism is completely general and it can be applied to other processes: the large logarithmic con- 
tributions are universal and, thus, they are systematically exponentiated in a process-independent 
form (see Eqs. ( TTUT) and (fT3l)); the process-dependent part is factorized in the hard-scattering coef- 
ficient ri. A constraint of perturbative unitarity is imposed on the resummed terms (see Eq. (TTT|) ) . 
so that the qx smearing produced by the resummation does not change the total production rate. 
This constraint reduces the effect of unjustified higher-order contributions at intermediate qr and 
facilitates the matching procedure with the complete fixed-order calculations at large qr- In par- 
ticular, when the rapidity is kept fixed, the integration over q T of da y '(dqTdy) at NNLL+NLO 
accuracy returns da/dy at NNLO. 

We have presented numerical results for Higgs boson production at the LHC. Comparing fixed- 
order and resummed calculations, we find that the resummation effects are large at small qr (as 
expected) and still sizeable at intermediate qT- The inclusion of the rapidity dependence has 
little quantitative impact on this picture since, as we have shown, the qT resummation effects are 
mildly dependent on the rapidity. Going from the central to the (moderately) off-central rapidity 
region, the qT shape of the spectrum slightly softens. In the range from small to intermediate 
values of q?, the residual perturbative uncertainty of the NNLL+NLO predictions for da/ (dqTdy) 
is comparable to that of advanced (NNLO or NNLL+NNLO) calculations of the qT inclusive cross 
sections da/dy and a. 



A Appendix 

In this appendix we present the structure of the resummation formula (flOj) by explicitly including 
the dependence on the flavour indeces of the colliding partons. 

In the context of our resummation formalism, a detailed derivation of exponentiation in the 
multiflavour case is illustrated in Appendix A of Ref. [33]. Considering a generic LO partonic 
subprocess c + c^F(F = H and c = c = g in the specific case of Higgs boson production by 
gluon fusion), and performing qT resummation after integration over the rapidity, the resummed 
component d&a^ /dq\ of the partonic cross section is controlled by the iV-moments W^ a2 N . The 
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final exponentiated result for these iV-moments is given by the master formulae (106)-(108) of 
Ref. [44]. We recall the master formula (106) in the following form: 

W a F ia2jJV (fo,M;a s ) = £ W^f* (M, a s ) exp{Q {IhN (a s , L)} , (18) 

{i} 

where the sum extends over the following set of flavour indices: 

{/} = c,c,i i ,i 2 ,b 1 ,b 2 ■ (19) 

and, for simplicity, the functional dependence on various scales (such as the renormalization and 
factorization scales) is understood. The functions G{i},n and Ti.^ F N are given in the master 
formulae (107) and (108), respectively 

In the present paper, qx resummation is performed at fixed values of the rapidity, and the 
double (Nx,N 2 ) -moments W^af 2)F in Eq. (ED replace the iV-moments Wf ia2jAr of Ref. |44J. The 
generalization of Eq. (TTUT) to the multiflavour case is straightforwardly obtained from Eq. ( fT8l) by 
the simple replacement N — > (Ni,N 2 ): 

Hiy )F (^;a s )^H2f'* )f W«s) «p{^^(as,L)} ■ (20) 

{i} 

The exponent Q^' N2) of the process-independent form factor and the process-dependent hard 
factor Hill,^ 1 '^ F are 

Gfn' N ^ = Qc + Qix,Nx + {?d>i,7Vi + Gi 2 ,N 2 + Qcb 2 ,N 2 i (21) 

'Hi I ? n '} Nl ' N2 " >F = o% j? S c C c bi.Ni -Ejvi Urn Ccbo.Ni Ek„ V M l U]\ 



a\a 2 



'N 2 



- 6 • (22) 

The expressions in Eqs. (1211) and fl22l) are completely analogous to the master formulae (107) and 
(108) in Ref. [44] (the functional dependence on the scales M, [Ir, \ip and Q is explicitly denoted in 
those formulae). In particular, we note that the dependence of C}( Nl ' N2 *) and 7-t( Nl ' N2 '> on the Mellin 
variables Ni and N 2 is completely factorized: each of terms on the right-hand side of Eqs. (|21"1) and 
(|2"2"|) depends only on one Mellin variable (either Ni or N 2 ). This factorized structure is completely 
consistent with Eq. f|T4]) and with the physical picture discussed below Eq. (fill) : the dependence 
on N\ (N 2 ) follows the longitudinal-momentum flow and the flavour flow a\ — > b\ — > %\ — > c 
(a 2 — > b 2 — > ^2 - c) that are produced by collinear radiation from the initial-state parton with 
momentum pi (p 2 )- The various Mellin functions (Gi,N, E^j ,Un and so forth) in Eqs. (1211) and 
(1221) can be found in Ref. H. 
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